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Spin susceptibility of Anderson impurities is a key quantity in understanding the physics of Kondo 
screening. Traditional numerical renormalization group (NRG) calculation of the impurity contri¬ 
bution Ximp to susceptibility, defined originally by Wilson in a fiat wide band, has been generalized 
before to structured conduction bands. The results brought about non-Fermi-liquid and diamagnetic 
Kondo behaviors in Ximp, even when the bands are not gapped at the Fermi energy. Here, we use 
the full density-matrix (FDM) NRG to present high-quality data for the local susceptibility xioc and 
to compare them with Ximp obtained by the traditional NRG. Our results indicate that those exotic 
behaviors observed in Ximp are unphysical. Instead, the low-energy excitations of the impurity in 
arbitrary bands only without gap at the Fermi energy are still a Fermi liquid and paramagnetic. We 
also demonstrate that unlike the traditional NRG yielding xioc less accurate than Ximp, the FDM 
method allows a high-precision dynamical calculation of xioc at much reduced computational cost, 
with an accuracy at least one order higher than Ximp. Moreover, artifacts in the FDM algorithm 
to Ximp, and origins of the spurious non-Fermi-liquid and diamagnetic features are clarified. Our 
work provides an efficient high-precision algorithm to calculate the spin susceptibility of impurity 
for arbitrary structured bands, while negating the applicability of Wilson’s definition to such cases. 

PACS numbers: 71.27.-|-a, 75.20.Hr, 05.10.Cc 


I. INTRODUCTION 

The understanding of quantum impurity systems is of 
prime importance in condensed matter physics. Such 
systems, consisting of a small subsystem (the impurity) 
with only few degrees of freedom coupled to a continu¬ 
ous bath of noninteracting particles, often exhibit rather 
complex behavior due to the nontrivial interplay between 
Coulomb repulsion at the impurity site and the impurity- 
bath interaction. A typical example is the Kondo effect 
[H describing the screening of a localized magnetic mo¬ 
ment by conduction electrons. Kondo screening mecha¬ 
nisms are mostly characterized by the spin susceptibil¬ 
ity of the impurity, a fundamental quantity describing 
the impurity’s magnetic response. Quantum impurities, 
such as magnetic molecules, qubits, adatoms, or quan¬ 
tum dots, would also constitute an enormous miniatur¬ 
ization of data processing technology by encoding and 
storing information in their magnetic states Q. Making 
this prospect feasible requires very accurate knowledge 
on the local magnetic response of the impurity. Exper¬ 
imentally, such local magnetic properties are relevant in 
traditional nuclear magnetic resonance and neutron scat¬ 
tering experiments [l[. 

Theoretically, a prototype model of quantum impuri¬ 
ties is the single-impurity Anderson model Q- While 
exact solutions of the Anderson model are available by 
the Bethe ansatz in some limiting cases, Wilson’s 


numerical renormalization group (NRG) provides 

a systematic nonperturbative method for arbitrary im¬ 
purity systems. For the impurity magnetic response in 
the Kondo physics, the NRG calculation of the Ander¬ 
son model often emphasizes Wilson’s definition Q of the 
impurity contribution Ximp to the susceptibility due to 
a global magnetic field (in the following, we term it the 
impurity susceptibility as usual), rather than the stan¬ 
dard definition @ of the local spin susceptibility xioc in 
response to a local field. Despite the experimental rele¬ 
vance, this is because evaluation of the local susceptibility 
is equivalent to calculating a dynamical spin correlation 
function which is much more complex and less accurate in 
the traditional NRG [m or its reduced density-matrix 
extension [^. Alternatively, one may calculate xioc by 
differentiating the local magnetization with respect to 
the external field S' From a numerical point of 
view, however, performing differentiations is something 
to avoid if possible. Since in the most relevant case of 
a flat wide band the impurity susceptibility Ximp exactly 
coincides with the local susceptibility xioc it is clearly 
more convenient to only calculate Ximp which is a ther¬ 
modynamic quantity and thus can be very accurately ob¬ 
tained using the traditional NRG. 

However, the nature of the electron bath in the Kondo 
problem varies from one realization to another. It may 
possess a band of strong energy dependence, as in mag¬ 
netic impurities adsorbed on graphene [HI and square 
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lattice [12|, doub le q uantum dots [l^, the narrow-band 
Anderson model 1^, and so on, for which the two sus¬ 


ceptibilities may differ significantly 1^. In this case, 


extreme caution should be exercised when applying Wil¬ 
son’s definition, because the impurity susceptibility, be¬ 
ing defined through global quantities, involves a subtlety 
from the conduction bath such that it may not necessarily 
yield the intrinsic spin dynamics of the impurity. Partic¬ 
ularly, a negative Ximp was found by the calculations 

signaling a diamagnetic Kondo impurity, even with 
non-Fermi-liquid behavior violating the low-temperature 
plateau in Ximp EH- This is inconsistent with one’s 
intuition. We emphasize that in those systems [liHii, 
though strongly energy dependent, the conduction bands 
are not gapped at the Fermi energy. Therefore, the im¬ 
purity Kondo physics is very likely to be still conven¬ 
tional, without the exotic non-Fermi-liquid and diamag¬ 
netic behaviors (gapped bands [l^ could, of course, lead 
to non-Fermi-liquid behavior). With this respect, in or¬ 
der to provide a definitive answer to whether an energy- 
dependent band without gap could indeed make the im¬ 
purity diamagnetic and a non-Fermi liquid, it is highly 
desirable to establish a dynamical method of precisely 
calculating the local susceptibility. This is because xioc is 
the standard quantity @ that characterizes the intrinsic 
spin dynamics of Anderson impurities in arbitrary con¬ 
duction bands. 

Recent NRG improvements of the full density-matrix 
(FDM) generalization Il7| based on the complete basis 
set of discarded states [18| may fulfill this purpose. This 
FDM approach avoids the overcounting ambiguity and 
the single-shell approximation for the density matrix in 
the previous NRG, giving rise to very accurate dynami¬ 
cal correlation functions for local operators at arbitrary 
temperature. Within the FDM algorithm, the zeroth- 
mome nt sp ectral sum rule holds exactly to machine pre¬ 
cision [H El- Higher-moment spectral sum rules are 
also fulfilled with a high accuracy [l^. This motivates 
us to wonder whether the local susceptibility, as a (—l)th 
moment of the impurity-spin correlation function, can be 
calculated accurately from the FDM method, with an ac¬ 
curacy even higher than the impurity susceptibility. Note 
that the accuracy of the impurity susceptibility from the 
FDM, as indicated by Ref. [10], is only in the percent 
range. Besides, at low temperature, there also exists a 
severe artifact in the FDM algorithm to Ximp (detailed 
in Sec. Ill A). 

In this paper, we apply the FDM approach to dynam¬ 
ically calculate the local susceptibility of an Anderson 
impurity, without involving numerical differentiations. It 
is shown that within the FDM algorithm, the local sus¬ 
ceptibility xioc can be calculated very accurately in all 
parameter regimes, while the thermodynamic calculation 
of the impurity susceptibility Ximp encounters severe ar¬ 
tifacts at low temperature. By comparison with exact 
Bethe ansatz results, the FDM results of the local sus¬ 


ceptibility Xioc are found accurate in the permil range. 
This accuracy is obtained at much reduced computa¬ 
tional cost, and is at least one order higher than the 
accuracy of the impurity susceptibility. We then revisit 
the previous studies [lll - ll3| . by applying the FDM algo¬ 
rithm for the local susceptibility xioc to the case with ar¬ 
bitrary energy-dependent conduction bands. The revisit 
reaches a definitive answer: as long as the host density 
of states is not gapped at the Fermi energy, the low- 
energy excitations of an Anderson impurity in arbitrary 
energy-dependent bands are still a Fermi liquid and para¬ 
magnetic. The resultant Kondo screening of the impu¬ 
rity magnetic moment is qualitatively same as in a flat 
wide band. For such systems, the exotic behaviors previ¬ 
ously observed [lll - ll4| in the impurity susceptibility do 
not represent the correct Kondo physics of the impurity. 
In particular, we demonstrate that the non-Fermi-liquid 
property and diamagnetism are spurious behaviors aris¬ 
ing from the additional susceptibility of conduction elec¬ 
trons, which may be vulnerable to the NRG discretization 
error. 


II. MODEL HAMILTONIAN AND FDM 
APPROACH TO SPIN SUSCEPTIBILITIES 

We consider the single-impurity Anderson model de¬ 
scribed by the Hamiltonian H = Hbath + + H-mt, 


.Hbath — 

k,(T 

-^imp — ^ ^ 

<7 

Hint = J2^kCi,d, + R.c.. 


( 1 ) 

( 2 ) 

( 3 ) 


where Cl^ (rft) creates an electron with energy et (£d) 
and spin cr =t,i in the bath (impurity), rio- = 

U parametrizes the on-site Coulomb repulsion, and the 
two subsystems are coupled via the hybridization 14. 
The influence of the bath on the impurity’s dynam¬ 
ics is fully determined by the hybridization function 
r(e) = |14p(5(e — Sk), which depends on specific 

realizations of the electron bath. In the standard case, 
F(e) is constant for \e\ < D and zero otherwise, with 
D being the half bandwidth. While D is usually the 
largest energy scale of the problem, the effect of finite 
bandwidth becomes important in the narrow-band model 
0. For a magnetic impurity adsorbed on the top of a 
carbon atom in graphene with Rashba spin-orbit inter¬ 
action, the resulting hybridization function of the model 
has a linear energy dependence with sharp discontinu¬ 
ities llj. r(e) may have even singularities around the 
Fermi energy, such as for magnetic impurities in square 
lattice [l^. When the single-impurity Anderson model 
pertains to the double quantum-dot system consisting of 
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an interacting dot coupled to the leads through another 
noninteracting dot, the impurity is coupled to an effec¬ 
tive bath with a Lorentzian hybridization function [l^ . 
Following, we calculate the local spin susceptibility yioc 
of the above Anderson model with arbitrary forms of r(£:) 
by using the FDM-NRG method. 

The NRG strategy starts from discretizing the 

bath spectrum on a logarithmic grid of energies ±T)A“" 
with A > 1 and n = 0,1, 2, • • •, thereby transforming the 
original impurity model into a semi-infinite tight-binding 
chain with exponentially decreasing hopping matrix ele¬ 
ments, via a standard tridiagonalization procedure after 
dropping high-mode states in each discretization inter¬ 
val. With its first site representing the impurity, the 
chain is then diagonalized iteratively adding one site at 
each step. In order to restrict the exponentially grow¬ 
ing Hilbert space, eigenstates of the chain Hamiltonian 
Hm including the newly added site (the A^th site), are 
constructed from the states of the A^th site and the Mk 
lowest-lying eigenstates (kept states) of the chain Hm-i 
without the A^th site, while discarding the remaining 
eigenstates of The iterative diagonalization pro¬ 

ceeds until the hopping matrix element between the last 
added site, say N = A^maxj and its immediate neighbor 
becomes the smallest energy scale of the problem, such 
that the Hamiltonian of the full chain represents 

a good approximation of the original Anderson model. 

Anders and Schiller [l^ have introduced a complete 
basis set of the Fock space of by constructing 

the tensor-product state \l,e-,N) = \l',N) 0 |aAr+i) 0 
law- 1 - 2 ) • • laWmax) from |/; A^) the 1th discarded states of 
Hpj and \am) the state of the mth site with am = {o,t4 
: ti}; where e denotes collectively the degrees of freedom 
of the sites m = N + 1, - ■ ■ , A^max, i-O., the environment 
of i?w- Let A^min being the first iteration at which high- 
energy states are discarded and taking all eigenstates of 
the last iteration A^max as discarded, the completeness 
relation of the basis set {|1, e; A^)} reads as 

^max 

^ ^|l,e;fV)(l,e;fV| = l, (4) 

iV=-/Vmin 

along with the orthonormality 

{l,e-,N\l\e'-N') =5iv5ee'5NN', (5) 

and a useful identity for the subspaces spanned by the 
kept states (denoted with k) at iteration N, and by all 
discarded states at succedent iterations N' > N: 

^max 

^|fc,e;Ar)(fc,e,iV|= ^ ^ |1, e; e; iV'|. 

k,e N'—N-\-l l,e 

( 6 ) 

Note also that for N' < N, the kept and discarded states 
are orthogonal (fc, e; A^|l, e'; A^') = 0. Since |s, e;A^) 
(s = 1, fc) is only an exact eigenstate of iJw correspond¬ 
ing to an eigenvalue with 4^“““^-fold degeneracy. 


one has to assume it is also an eigenstate of the original 
model, H\s,e-,N) Ri E^\s,e-,N). This so-called NRG 
approximation represents the only approximation of the 
FDM algorithm. Weichselbaum and Delft [13 hence 
write the full density matrix p of iL as follows: 

1 ^max 

P=Z^ ^e-''<|/,e;fV)(/,e;fV|, (7) 

N—Nmin l,e 

^max 

W=Wmi„ I 


with P = 1/{ksT) the inverse temperature. By using this 
form of the density matrix, the complete basis set of dis¬ 
carded states, and the NRG approximation, all dynamic 
and static properties of El can be evaluated. 

By definition, the local susceptibility describes the im¬ 
purity magnetization in response to a weak magnetic field 
B applying only at the impurity site Q 

_ d{Mm)H+H'^ I \2 A 1C I \C \ 

Xioc = hm -- = spb) / dr {S,{t)S,}h, 

B^o dB Jq 

( 9 ) 

where Mm = gfJ-B^z is the impurity magnetization op¬ 
erator, g the Lande g factor, g,B the Bohr magneton, 
Sz = — is the 0 component of impurity spin, and 

S'j;(t) = Sze~'^^. {■ ■ ■ )h+h^ and ■ ■)h denote the 

thermodynamic average with respect to the Hamiltonian 
E[ with and without the perturbation El'm = —ggis^zB. 
The second equality in Eq. (9), used {Sz)h = 0, rep¬ 
resents an exact mathematical relation which expresses 
a response in the impurity magnetization due to an in¬ 
finitesimal local field in terms of an imaginary-time Mat- 
subara Green’s function [^. This already constitutes an 
operational ground for dynamical calculating the local 
susceptibility, while obviating the need to evaluate a nu¬ 
merical derivative. However, the FDM algorithm for the 
Matsubara function would involve the NRG approxima¬ 
tion, e; A^) ~ |s,e;N), which at low tem¬ 

perature becomes severe [i^ because of large r involved. 

To obtain high-quality susceptibility data in the low- 
temperature Kondo regime, it is thus better to work with 
the retarded Green’s function for which the quality of the 
NRG approximation does not rely on the temperature. 
Using the Kubo formula for linear response in the static 
limit, the local susceptibility can be rewritten as 


Xloc 


Gs, (e) 


-{gfiB)'^Re[Gs^{0)] 

TT J -00 ^ 

- dte7^^‘0(t)([54t),At(O)])w. 


( 10 ) 

( 11 ) 


It is easy to verify the equivalence of Eqs. (9) and (10) 
from their Lehmann representations. Moreover, the sec¬ 
ond equality in Eq. (10) reveals that the local suscepti¬ 
bility is just the (—l)th moment of the impurity-spin 
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spectral function. To evaluate within the FDM algo¬ 
rithm 0 the retarded Green’s function of the local spin, 
Gs^{e), the full density matrix Eq. (7) and the complete 
basis set Eq. (4) are inserted into the thermal average 
in Eq. (11). Making use of the properties of the basis 
[Eqs. (5) and (6)] and taking the NRG approximation 
e; iV) Ri ‘|s, e;A^), we end up with (see 

Appendix for details) 


Xloc 


_ 'ST' nG) r,NN' rc ]N rc ]N 

— 2^ '^Nkl ' Pk'k ■ Pzjfei ■ [^z\lk' 

AT, —A^min k,k' 

A/^max — 1 

+ E 

A/^ —A/^min k,l 
max 

+ E '^^NW ' l^2]u' ' (12) 

JV=Ar„i„ 1,1' 


Here restricts the summation to N < N', the matrix 
elements of local spin = (s; A^), and the 

reduced density matrix = Yei^'^\PN'\k,e-, N) 

in which the density matrix of shell N', 


PN' = (13) 

l,e 


satisfies p = Pn- Finally, the three coefficients 

in Eq. (12) are given by 


used quantity is defined originally by Wilson [^, to be 
the difference of the total magnetic response with and 
without the impurity, under a global field B, 


Ximp 


d 

hm — 
B-s-o oB 


+ (Afb)B+B' 


-(A7b)Bbath-rff^ 

= {gpB?P{{St)H-{Sl)H.^Y, (17) 


where H'^ = -gpuStB, Mb = gpsSb, = -gpuShB, 
and St = Sb + Sz^ with Sb being the z component 
of bath spin. The second equality, using [S't, H] = 
[S';,, Mbath] = 0, indicates that Ximp is a thermodynamic 
quantity which, as usual, can be obtained accurately us¬ 
ing the traditional NRG. Within the FDM method, at the 
A'th iteration, both St and Sb should be divided into (i) 
the z-component spin, Sn■, of the resultant chain Hamil¬ 
tonian Hjsi with or without the impurity site, and (ii) 
the z-component spin of iVniax ~ N environmental sites. 
Following Ref. [10], the FDM formula for Ximp is 


Ximp 

PigpB) 

-4 


_ _ ^with impurity _ jj^without impurity 

A/^max 

W=Afmin I 


+ (A^max 



(18) 


(19) 


— Fp 

'-'Nkl — 

^(2) _ _ 

'-'Nkl ~ 


-E^)+ip^ 


1 A^max — AT 


^(3) _ 

^NW ~ 


e-^Erc(^) 


^ - ~ NkE 

4A7^ax-A7 

^ P{Ef-EP) ’ 


(14) 

(15) 

(16) 


where ry is a dimensionless infinitesimal (already absorbed 
/3) to deal with the accidental degeneracy of kept and 
discarded states, E^ = E^, which appears when high- 
energy states are truncated at a degenerate eigenenergy 
of Mat. Keeping p fixed here is equivalent to broaden 
discrete 6 functions in the spectral function lm[Gs^(w)] 
using a Lorentzian kernel with temperature-dependent 
width {p/j3). In this work, we always set p = 0.001, un- 
less stated otherwise. Since the value of in the limit 
E^ = E[f is well defined, there is no need to introduce an 
infinitesimal imaginary part to Eq. (16). After perform¬ 
ing a “forward run” along the Wilson chain to iteratively 
generate all relevant NRG eigenenergies E^, eigenstates 
|s; A^), and matrix elements we can evaluate the 

reduced density matrix in a single “backward run” 
[13 , thereby obtaining all information needed for dynam¬ 
ically calculating the local susceptibility. 

For comparison, we also present the FDM results for 
the impurity contribution to susceptibility. This widely 


At first glance, the absence of the reduced density matrix 
in Eq. (19) implies that thermodynamically evaluating 
Ximp by two NRG runs (one with and one without the im¬ 
purity) would require less computational resources than 
dynamically evaluating xioc by a “forward” and “back¬ 
ward” NRG run. As we show in the following, however, 
within the FDM approach the accuracy obtainable for 
Xloc even at lower computational cost can be one order 
higher than that of Ximp- 

To establish a relation between the local xioc and impu¬ 
rity Ximp susceptibilities, let us consider their difference 
Xloc — Ximp- From the equation of motion (EOM) of the 
bath-electron Green’s function, one readily obtains 


Xloc-Ximp = 5PB|im — 


ds 
27r 

,d^aie) 


a 

-GdAe) 


ds 


6G{e) 

( 20 ) 


where f{e) stands for the Fermi-Dirac function, ^^(e) = 
Ek |t4P/(e- Efc + ^agpsB + i0'^) is the self energy due 
to the impurity-bath coupling, and SG{e) = Gd„{s) — 
G'^ (e), with Gd,„ (e) [G'^ (e)] being the impurity retarded 
Green’s function corresponding to the Hamiltonian H 
in the global (local) magnetic field. The term contain¬ 
ing Gd„{s)dYj^{s)/de in Eq. (20) represents the addi¬ 
tional susceptibility from the conduction electrons M- 
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For a flat band in the wide-band limit, one has exactly 
5So-{e)/9£ = 0 and SG{e) = 0, leading to yioc = Ximp 
as expected from the Clogston-Anderson compensation 
theorem [^. On the other hand, the breakdown of this 
theorem in the presence of a narrow bandwidth and/or 
strong energy dependence in So-(£) [or r(e), equivalently] 
0 can result in yioc and Ximp differing substantially. In 
this case, the impurity susceptibility may be no longer 
suitable for faithfully characterizing the intrinsic spin dy¬ 
namics of the impurity, especially when the additional 
bath susceptibility becomes significant. 


III. RESULTS AND DISCUSSIONS 

What follows are the numerical results calculated by 
the FDM NRG in the units of D = gfiB = fcs = 1- Tech¬ 
nically, the NRG is still an approximation method that 
involves the discretization error (controlled by the dis¬ 
cretization parameter A) and the truncation error (con¬ 
trolled by the number of kept states in each iteration 
and also A) [2^ . The two types of errors are interrelated. 
For coarser discretization at larger A, the discretization 
error increases, whereas the truncation error decreases 
due to the enhanced separation of energy scales. Large 
A may also introduce spurious oscillations into thermo¬ 
dynamic quantities. These oscillations can be removed 
by using the z-averaging procedure, wherein one averages 
the final results from independent NRG calculations for 
Nz interleaved discretization meshes 
with Nz values of the twist parameter z equally dis¬ 
tributed in [0,1). As for Mk, unlike in the conventional 
method, the number of kept states needed by the FDM 
algorithm can be largely reduced while still obtaining sat¬ 
isfactory accuracy for physical observables, due to the use 
of a complete basis set. Particularly, the sum-rule nature 
inherent in the FDM algorithm to yioc makes the local 
susceptibility very insensitive to Mk- We thus always 
reduce the number of kept states for calculating yioc to 
half the number of states kept for Ximp- Nevertheless, 
one still has to carefully choose these parameters in any 
practical NRG calculations such that the resultant dis¬ 
cretization and truncation errors do not affect physical 
conclusions drawn from the NRG data, which should be 
robust to changes of parameters. 


though the two susceptibilities are essentially identical 
in this case, one can not expect the FDM-NRG calcu¬ 
lation would always yield equal results for them. This 
is because the errors inherent in NRG may affect xioc 
and Ximp in different ways since different algorithms as 
Eqs. (12) and (18) are adopted. For small A and moder¬ 
ate Mk [Fig. 1(a)], it is shown that at high temperatures 
both algorithms have introduced only tiny artifacts. But 
a severe artifact appears in the low-temperature behavior 
of Ximp, which obviously breaks the Fermi-liquid prop¬ 
erty. At low temperatures, the Fermi-liquid theory [l| of 
a Kondo impurity demonstrates a linear temperature de¬ 
pendence of the effective Gurie constant Ty, i.e., the spin 
susceptibility x would develop a plateau as T —)> 0. This 
plateau is indeed very well formed in yioc as a function 
of T, indicating that dynamical FDM calculation of the 
local susceptibility is much more accurate than the ther¬ 
modynamic calculation of the impurity susceptibility. 



FIG. 1: (Color online) Local xioc and impurity Ximp suscepti¬ 
bilities calculated by the FDM NRG for the Anderson model 
with a flat wide band, over a wide temperature range con¬ 
tinuously evolving from the high-temperature free-orbital to 
the low-temperature Fermi-liquid regimes. Model parame¬ 
ters: D — 1, r = 0.001, U/T = 12, and Sd/F = —6, unless 
indicated otherwise. Note that the curves in (c) regarding 
different parameters are highly coincident. NRG parame¬ 
ters in (d): for xioc (lines), A = 4, Mk = 300; for Ximp 
(symbols), A = 8, Mk = 600. The z averaging is per¬ 
formed only for Ximp in (b) and (d), using Nz = 8, z = 
0,0.125, 0.25, 0.375, 0.5, 0.625, 0.75, 0.875. 


A. Impurities in the flat wide conduction band 

We first examine the FDM-NRG calculations of the 
two susceptibilities for a flat conduction band [r(e) = F] 
in the wide-band limit [D F, 1/, jedj). The calcula¬ 
tions are performed for the temperature dependence of 
the susceptibilities, by fixing F while varying the on¬ 
site Coulomb repulsion and the impurity level position, 
as shown in Figs. 1 and 2 and Tables I and II. Even 


The artifacts uncovered in Fig. 1(a) are mostly due to 
the truncation errors which can be largely reduced by in¬ 
creasing the discretization parameter A. For the impurity 
susceptibility Ximp, we find that using a discretization 
parameter as large as A = 8 and performing the z aver¬ 
aging for Nz = 8 meshes are necessary in order to obtain 
satisfactory data [Fig. 1(b)], in agreement with Ref. [10]. 
On the other hand, high-quality data of xioc robust to 
changes of the number of kept states Mk and the param- 
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eter rj, are already available at moderate A = 4, as given 
in Fig. 1(c). Since at this value of A the local suscepti¬ 
bility xioc exhibits no spurious oscillations, there is no 
need to carry out the z averaging. This significantly low¬ 
ers the computational cost. With these NRG parameters 
tuned independently for yioc and Ximp to eliminate the 
artifacts, the two susceptibilities plotted as Ty vs T in 
Fig. 1(d) seem indeed identical in all parameter regimes. 
The resultant curves exhibit features of typical Kondo 
screening for —e^, Sd + U F, that is, Ty first increasing 
from the high-temperature value 1/8 of a free impurity 
towards its local-moment value 1/4 for intermediate tem¬ 
perature, and then falling to zero as T —^ 0 due to the 
screening of the local magnetic moment by conduction 
electrons. 

However, a severe drawback of the FDM method to 
Ximp, obscured by the very small values of Tximp at 
low temperature and thus overlooked by Ref. [10], is that 
the low-temperature artifacts appearing in yimp can not 
be completely eliminated in all parameter regimes, even 
though large A and the z averaging are used. To ex¬ 
plicitly demonstrate this, we present by Fig. 2 a detailed 
comparison of the two susceptibilities, plotted as y vs T 
rather than Ty vs T, emphasizing the low-temperature 
Fermi-liquid plateau. As depicted, our FDM calculation 
of the local susceptibility gives high-quality data for all 
impurity parameters, ranging from the strongly corre¬ 
lated to the noninteracting regime [Fig. 2(a)], and from 
the Kondo to the mixed valence and into the empty- 
orbital regime [Fig. 2(b)], whereas the data of the impu¬ 
rity susceptibility become wildly irregular at extremely 
low temperature. We have checked that such kind of ir¬ 
regularities always show up in yimp by further increasing 
A, A^ 2 , and even Mk- Seemingly, these irregularities are 
small in the Kondo regime and in the particle-hole sym¬ 
metry case, but becomes severe in the mixed-valence and 
empty-orbital regimes. The precise origin of these irreg¬ 
ularities is not clear at present. It may be 23| due to an 
imperfect cancellation of the logarithmic discretization 
oscillations by the simple z-averaging procedure, or the 
limitations in numerical precision for calculating yimp by 
taking the subtraction of two extensive macroscopic val¬ 
ues to obtain an impurity-related finite quantity. In any 
case, the FDM gives indisputably better data quality for 
the local susceptibility than for the impurity susceptibil¬ 
ity, being distinct from the conventional NRG. 

To quantitatively demonstrate the advantage of our 
FDM algorithm for the local susceptibility, we compare 
both susceptibilities with the exact zero-temperature re¬ 
sults xba from the Bethe ansatz Q , as shown in Tables I 
and II. Note that one can not take exactly T = 0 in FDM- 
NRG calculations. The zero-temperature results within 
NRG are actually extracted from the corresponding low- 
temperature results which must be convergent as T —>■ 0. 
This poses no difficulties for the local susceptibility, be¬ 
cause xioc do indeed converge to a definite value as T —>• 0 





lines: Xioc 
symbols: x,„p 

2Ej+U = 0 

u/r=i2(top), n, 

10 , 9 , 8 , 7 , 
6 , 5 , 4 , 3 , 2 , 
l,0(below) 


lO''^ 10'" 10'' 10'' lO'^' 10- 



FIG. 2: (Color online) FDM-NRG results of the local and 
impurity susceptibilities plotted as x vs T, focusing on the 
low-temperature Fermi-liquid plateau, for the symmetric (a) 
and asymmetric (b) Anderson models with a flat wide band 
{D = 1, F = 0.001). NRG and z-averaging parameters are 
the same as in Fig. 1(d). 


(see Fig. 2). We thus take the zero-temperature value of 
Xioc at the lowest temperature involved in our calcula¬ 
tion, i.e., T = 10“®. However, due to its low-temperature 
artifacts, extracting the zero-temperature value of ximp 
would be problematic since Ximp(T) is not convergent as 
T —>■ 0. To avoid the problem, the zero-temperature im¬ 
purity susceptibility we used in Tables I and H for com¬ 
parison is actually the values of Ximp(T') in the plateau 
region at T = 10“^To, where Tq is the low-energy Kondo 
scale defined by ToXimpiTo) = 0.0701 Comparison 
with the Bethe ansatz results indicates a relative error of 
Xioc from the FDM in the permil range, while Ximp is ac¬ 
curate only to within a few percent. This much improved 
accuracy of xioc becomes even more remarkable when 
considering the following fact. Our dynamical calcula¬ 
tion of the local susceptibility does not need to perform 
the z-averaging, and includes only half of the eigenstates 
as kept in thermodynamically evaluating Ximp, and thus 
is carried out at much reduced computational cost. 


B. Impurities in energy-dependent bands 

We now turn to clarify the effect of energy-dependent 
conduction bands on the magnetic response of an An¬ 
derson impurity. The single-impurity Anderson model 
with a structured conduction band has already been in- 
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TABLE I: Zero-temperature local xioc and impurity Ximp sus¬ 
ceptibilities from the FDM NRG in comparison with the ex¬ 
act result xba from the Bethe ansatz, for the symmetric An¬ 
derson model with a fiat wide band [D — l(NRG), oo(BA), 
r = 0.001, Ed = —U/2]. NRG and z-averaging parameters are 
the same as in Fig. 1(d). 


u/r 



% error 

^imp 

% error 


(smb)^ 

(smb)^ 

12 

9970.8 

9974.0 

-tO.03% 

10136 

+ 1.7% 

11 

6950.4 

6943.9 

-0.09% 

7063.7 

+ 1.6% 

10 

4854.3 

4854.8 

-tO.01% 

4931.4 

+1.6% 

9 

3397.7 

3398.3 

-tO.02% 

3451.3 

+1.6% 

8 

2383.8 

2382.7 

-0.05% 

2420.7 

+1.5% 

7 

1676.9 

1673.9 

-0.18% 

1702.5 

+1.5% 

6 

1183.2 

1182.3 

-0.08% 

1201.6 

+1.6% 

5 

837.65 

837.13 

-0.06% 

849.94 

+1.5% 

4 

595.23 

594.40 

-0.14% 

603.39 

+1.4% 

3 

424.76 

423.41 

-0.32% 

430.15 

+1.3% 

2 

304.54 

304.01 

-0.17% 

308.14 

+1.2% 

1 

219.50 

219.68 

-tO.08% 

221.54 

+0.9% 

0 

159.15 

159.09 

-0.04% 

159.34 

+0.1% 


TABLE II: Zero-temperature local xioc and impurity Ximp 
susceptibilities from the FDM NRG in comparison with the 
exact result xba from the Bethe ansatz, for the asymmetric 
Anderson model with a flat wide band [D = l(NRG), oo(BA), 
r = 0.001, [//F = 12]. NRG and z-averaging parameters are 
the same as in Fig. 1(d). 


sd/r 

^BA 

^loc 

% error 

^imp 

% error 



(9Mb)^ 

-5 

8748.3 

8749.9 

+0.02% 

8811.6 

+0.7% 

-4 

5910.2 

5904.8 

-0.09% 

5974.5 

+ 1.1% 

-3 

3078.5 

3075.9 

-0.08% 

3108.2 

+ 1.0% 

-2 

1246.7 

1244.9 

-0.14% 

1285.5 

+3.1% 

-1 

412.21 

412.13 

-0.02% 

391.54 

-5.0% 

0 

131.75 

131.42 

-0.25% 

140.54 

+6.7% 

1 

50.045 

49.996 

-0.10% 

49.232 

-1.6% 

2 

23.607 

23.563 

-0.19% 

23.084 

-2.2% 

3 

13.150 

13.128 

-0.17% 

13.294 

+1.1% 

4 

8.2236 

8.2098 

-0.17% 

8.5423 

+3.9% 

5 

5.5769 

5.5774 

+0.01% 

5.8513 

+4.9% 


tensively studied in the literature. It has been shown 
that when there is a gap (hard or soft) at the Fermi en¬ 
ergy in the host density of states, the impurity ground 
state can undergo quantum phase transitions from the 
Kondo screening to the local moment state, giving rise 
to non-Fermi-liquid behavior [l6| . However, there is also 
a class of systems [lll - [l^ in which the host density of 
states, although strongly energy-dependent, has i io g ap 
at the Fermi energy. Previous investigations |lll4l4l | of 
such systems suggest that the low-energy excitations of 
the impurity can still be a non-Fermi liquid and even dia¬ 


magnetic. The observations are quite surprising and are 
based on the results of the impurity susceptibility Ximp 
calculated by the traditional NRG. Since Ximp can not 
directly reflect the magnetic states of the impurity, we 
thus revisit this problem by directly calculating the local 
susceptibility xioc of these systems using the FDM NRG, 
in order to clarify whether an energy-dependent band 
without gap could indeed render the impurity diamag¬ 
netic and a non-Fermi liquid, or the previous observation 
[nnii is just spurious. 

The first system ll| we shall examine is a magnetic 
impurity adsorbed in graphene with Rashba spin-orbit 
interaction, in which the interplay of the Rashba cou¬ 
pling and the linear graphene dispersion results in an 
effective host density of states described by the following 
hybridization function IT| 


F(e) = Fo[|£| + a + (|e| - A)0(|£| - 2X)]/D, (21) 


where A characterizes the magnitude of the Rashba in¬ 
teraction, the prefactor Fq = floDV'^ / {4vp) with Hq the 
graphene unit-cell area, V the overlap between the im¬ 
purity level and the nearest carbon pz orbital, and vp 
the Fermi velocity. This hybridization function has a lin¬ 
ear energy dependence with discontinuities at e = ±2A. 
When the Fermi energy p of this system is tuned to lie ex¬ 
actly at the discontinuity p = 2A, Ref. [11] has found by 
the traditional NRG that lirnr^o^Ximp is not zero but 
rather negative, i.e., limT^o Ximp is negatively divergent 
by 1/T. This implies a non-Fermi-liquid and diamagnetic 
behavior in the impurity ground state. For the Fermi en¬ 
ergy not exactly at (but very close to) the discontinuity, 
while the Fermi-liquid property limT^-o T Ximp = 0 is re¬ 
stored, the impurity still passes through a ternperature 
window of diamagnetic behavior (ximp < 0) [ll|. We 
have calculated by the FDM NRG the local and impu¬ 
rity susceptibilities of this system. Results are presented 
in Fig. 3. Note that our FDM results of Ximp [Figs. 3(c) 
and 3(d)] do indeed verify the results of Ref. [11], despite 
the fact that the FDM algorithm for Ximp introduces arti¬ 
facts at extremely low temperature. As explained already 
in Section III.A, these artifacts are invisible when plot¬ 
ted as Tximp vs T [Fig. 3(c)], but apparently show up in 
Ximp vs T [see the inset of Fig. 3(d), where the artifacts 
have violated the Fermi-liquid plateau at T < 10“^° in 
the curves corresponding to /r = 2A ± 10“®]. 

However, the above non-Fermi-liquid and diamagnetic 
behavior found in Ref. [11] are not supported by our 
FDM results of the local susceptibility xioc- Figures 
3(a) and 3(b) demonstrate that for different values of 
the Rashba parameter, the low-energy excitations of the 
impurity are always a Fermi liquid (limT^-o Fxioc = 0, 
limT^o Xioc = const.) and paramagnetic (xioc > 0), no 
matter the Fermi level lies exactly at (/r = 2A) or slightly 
deviates from (p, = 2A ± 10“®) the discontinuity. More¬ 
over, unlike the behavior of Ximpj there is no significant 
difference in xioc for these positions of the Fermi level 
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FIG. 3: (Color online) Local [(a),(b)] and impurity [(c),(d)] 
susceptibilities, plotted as Tx vs T [(a),(c)] and x vs T 
[(b),(d)], from the FDM for a magnetic impurity in graphene 
with different Rashba spin-orbit coupling A and chemical po¬ 
tential The resulting Anderson model 0 is characterized 
by the hybridization function Eq. (21). The legend in (a) is 
applied to all figures including the insets. Insets: impurity 
susceptibility for different jj, at fixed A = 0.006, plotted as 
Tximp vs T [inset of (c)] and yimp vs T [inset of (d)]. Model 
parameters: U = —2sd = O.SFq = 0.02, being the same as in 
Figs. 5 and 7 of Ref. [11]. NRG parameters: A = 2.5; for Xioc, 
Mk = 600, without the z averaging; for Ximp, Mk ~ 1200, 
Nz = 2 with 2 = 0,0.5. 


(see the red solid lines and the symbols in Fig. 3). This 
behavior of xioc is consistent with the underlying Kondo 
physics, while Ximp is not. Generally speaking, at tem¬ 
peratures much lower than the Kondo scale T Tq, 
only those conduction electrons within the energy win¬ 
dow |e — /r| < To participate in the Kondo screening. 
Since the Fermi energies ^ = 2A, 2A ± 10“® produce the 
almost same Kondo scale m which is far, far larger 
than the energy difference in these the same portion of 
conduction electrons around the Fermi level are involved 
in screening the impurity spin. Consequently, for these 
Fermi energies, the impurity magnetic response should 
also be almost equal, as indicated by our xioc rather 
than Ximp in Ref. [11] suggested. We thus argue that the 
correct magnetic property of the impurity in graphene 
with the exotic hybridization function Eq. (21) is still a 
standard Fermi liquid. 


The second system [l^ we have revisited is a magnetic 
impurity in the two-dimensional square lattice, with the 
half bandwidth D determined by the nearest-neighbor 
hopping energy D = it. Its host density of states p{e) 
has a Van Hove singularity near the Fermi energy fi = 0. 
According to Ref. [12], the distance A from the singu¬ 
larity to the Fermi energy is A = it' with t' the next 


nearest-neighbor hopping, and 


^ 21n[(4yj^^)/(|g + A|)] 

^ 7r2Dv'l- (A/D)2 

Based on the NRG results of Ximp within the Kondo 
model. Ref. [12] predicted a non-Fermi-liquid and dia¬ 
magnetic regime at low temperature for A = 0, in which 
RXimp ~ —0.072/1 ln(T/Zl)|°-^’^ (also leading to divergent 
limr^o Ximp)- For nonzero but very small A, the Fermi- 
liquid behavior Tyimp = cT is restored with the scale 
factor c remaining negative. These predictions are qual¬ 
itatively verified by our FDM results of Ximp based on 
the Anderson model, as shown in Fig.4(a). Our FDM 
results of xioc [also presented in Fig. 4(a)] again do not 
support these non-Fermi-liquid and diamagnetic behav¬ 
iors. While the authors of Ref. [12] attributed these spuri¬ 
ous behaviors to an overcompensation of the local spin by 
the conduction electrons, we draw a conclusion from the 
local susceptibility that a Kondo impurity in the square 
lattice with Van Hove singularities is still a Fermi liquid 
and there is no diamagnetic or overscreening effects. Ad¬ 
ditionally, for the values of A used in plotting Fig. 4(a), 
there are no sizable deviations in the local susceptibility. 
This is due to the same reason as already explained in 
the first system, i.e., these values are very close, which 
give rise to the nearly same Kondo temperature and thus 
involve the same portion of conduction electrons in the 
Kondo screening. 

We have also examined the effects of a Lorentzian and a 
narrow conduction band on the impurity spin dynamics. 
The Lorentzian host density of states can be realized by 
double quantum dot systems, where an interacting quan¬ 
tum dot (the impurity) is indirectly coupled to the leads 
through a noninteracting quantum dot 0. When the 
resonance in the noninteracting dot is at the Fermi en¬ 
ergy, the double quantum dot system maps onto a single¬ 
impurity Anderson model having a Lorentzian hybridiza¬ 
tion function 0 


r(£) 


Ai 

A2 £2 + A2 ’ 


(23) 


where A is the interdot coupling and A 2 the coupling 
between the noninteracting dot and the leads. On the 
other hand, the narrow-band Anderson model is gener¬ 
ated by the dynamical mean field theory for the Mott- 
Hubbard transition in infinite dimensions 1^ . Although 
the two models are of interest in such different contexts, 
the resulting impurity properties are similar since the 
Lorentzian band is in analogy with the narrow band in 
the sense that it diminishes the high-energy conduction 
states reducing effectively the bandwidth. Previous stud¬ 
ies 0,0] within the traditional NRG show that for both 
models the impurity susceptibility suggests a Fermi liq¬ 
uid (limT^oFXimp = 0, limT->oXimp = const.) in all 
parameter regimes, but also a diamagnetic (ximp < 0) 
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T 


FIG. 4: (Color online) Local and impurity susceptibilities 
from the FDM for (a) a magnetic impurity in square lat¬ 
tice [12], (b) the double quantum dot system [l3| . and 
(c) a narrow-band system M- The resulting hybridiza¬ 
tion function in the Anderson model is (a) r(£) = tvV'^ p{£) 
with p(e) given by Eq. (22), (b) given by Eq. (23), and (c) 
r(e) = r0(D — |£|) for D ^ ?7, |£d|. Model parameters: 
(a) U = —2ed = 0.05, 2V‘^ I{'kD) = 0.0006; being consis¬ 
tent with the Kondo-model parameters used by Ref. [12]; (b) 
U = —2£d = 0.5, A 2 = 0.02, being the same as in Fig. 2(d) of 
Ref. [13]; (c) U = —2ed, F = 10, being the same as in Fig. 1 
of Ref. [14]. NRG parameters: A = 2; for xioc, Mk ~ 600; for 
Ximp, Mk = 1200; without the z averaging. 


region in Ximp vs T for some parameters. As expected, 
our Ximp by the FDM repeats these properties [see the 
dashed curves in Figs.4(b) and 4(c)]. The Fermi-liquid 
property is also confirmed by our local susceptibility xioc 
presented in Figs. 4(b) and 4(c). However, xioc for both 
models is always positive even in the parameter regimes 
where Ximp is negative. This rules out the scenario of dia¬ 
magnetic impurities caused by the Lorentzian or narrow 
conduction bands. 


C. Origins of the spurious diamagnetic and 
non-Fermi-liquid behaviors in Ximp 


To pinpoint the origin of the diamagnetism and non- 
Fermi-liquid behavior in the impurity susceptibility, let 
us look back into the definition Eq. (17) of Ximp- The 
first term in Eq. (17) is always pos¬ 

itive and its contribution is very similar to the local 
susceptibility xioc- Their difference gives rise to the 
6G{e) term in Eq. (20). The other two terms in Eq. (17), 
\imB^o-^{{Mb) h+h^ -(Mh )represents 
the additional susceptibility of conduction electrons in¬ 
duced by the presence of the impurity. 6xc is exactly the 
Gd^ {e)dTicy{e)/de term in the right-hand side of Eq. (20). 
Due to the derivative of the self energy, the sign and 
magnitude of 5xc are very sensitive to the shape of the 
conduction band. It is this additional bath susceptibility 
included in the definition of Ximp that becomes negative 
and divergent at low temperature when the host den¬ 
sity of states is strongly energy dependent. Therefore, 
the diamagnetism and non-Fermi-liquid behavior found 
in previous studies [11-14] does not directly reflect the 
intrinsic impurity properties, and has nothing to do with 
the Kondo screening of the local moment at low temper¬ 
ature. 

For a deep insight into the non-Fermi-liquid property 
suggested by Ximp in the graphene [ll| and the square lat¬ 
tice 12| systems, we consider the U = 0 Anderson model 
which allows to clarify wether the Ximp ^ T dependence 
in the U ^ 0 case is qualitatively different from the ?7 = 0 
case and also allows to compare the FDM-NRG results 
with the exact ones. In the noninteracting case, the exact 
local and impurity susceptibilities can be obtained using 
the EOM approach, 


Aio 


Ximp 




(24) 




+ 


e- Ed - S(£) 

' i-is(£) 


1 2 


S(£) = 


de' 


e-Ed- E(£) 

r(e') 


E — e' + 10+ 


(25) 

(26) 


Note that for arbitrary energy-dependent bands un¬ 
gapped at the Fermi energy, Eq. (24) always give rise to 
finite (not divergent) values of xioc, as the temperature 
T —>• 0. But this is not the case for Eq. (25) of Ximp. 
Figure 5 presents the FDM-NRG results for the spin 
susceptibilities of the noninteracting Anderson impurity 
adsorbed in graphene [Fig. 5(a)] with the hybridization 
function Eq. (21) and in the square lattice [Fig. 5(b)] with 
the density of states Eq. (22). These are in good agree¬ 
ment with the exact EOM results. It is demonstrated by 
Fig. 5 that even in the U = 0 case, the impurity suscep- 
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FIG. 5: (Color online) Local and impurity susceptibilities 
from the FDM, for a noninteracting impurity (a) in graphene 
with the hybridization function Eq. (21) and (b) in the square 
lattice with the density of states Eq. (22). Model param¬ 
eters: (a) ?7 = £d = 0, Fo = 0.04; (b) U = Sd = 0, 
2 V‘^/{'!tD) = 0.0006. NRG parameters are the same as in 
Figs. 3 and 4, respectively. The exact EOM results, calcu¬ 
lated from Eqs. (24)-(26), are also presented for comparison. 


tibility Ximp is already negatively divergent as T —>-0, 
being qualitatively analogous to the corresponding in¬ 
teracting systems [see Figs. 3(c), 3(d), and 4(a)]. This 
qualitative analogy between the interacting and nonin¬ 
teracting systems confirms again that the graphene 11 j 
and square lattice [l^ systems are indeed a Fermi liquid 
[ij . From this point of view, the non-Fermi-liquid physics 
proposed previously 12| according to the divergence of 
Ximp as T —>■ 0 is conceptually incorrect and misleading. 

Aside from being vulnerable to produce misleading re¬ 
sults due to the involved Sxc term, the NRG calculation 
of the impurity susceptibility Ximp has another poten¬ 
tial flaw in a more fundamental aspect. Note that the 
NRG logarithmic discretization is an approximate proce¬ 
dure transforming the continuous conduction bath into 
the discretized Wilson chain. The bath properties are 
qualitatively changed by this procedure due to the dis¬ 
cretization error. With this respect, the NRG method is 
only suitable to calculate local quantities which do not 
explicitly involve the bath degrees of freedom, e.g., xioc- 
The NRG calculation of any nonlocal quantity explicitly 
involving the bath degrees of freedom, e.g., Sxc and thus 
Ximp, may be not reliable. A representative example is 
the spin susceptibility Xc = (5Ms)^/3('S'h)i/bath of conduc¬ 
tion electrons [i.e., the last term in Eq. (17)] in a flat wide 
band [p(£) = Pg], as shown in Fig. 6. For the original 



FIG. 6: (Golor online) Spin susceptibility Xc of conduction 
electrons in a fiat wide band. For the original continuous 
model, Xc is the temperature-independent Pauli susceptibil¬ 
ity. For the discretized Wilson chain, Xc is evaluated using 
the exact diagonalization method, and presented for Xc/Nt 
with Nt the total number of sites in the chain. We choose a 
sufficiently large number of sites to ensure that the coupling 
between the last two sites, ~ is far less than the 

temperature. 


continuous model, Xc gives the temperature-independent 
Pauli paramagnetic susceptibility Xc = Po- But 

the corresponding quantity in the discretized Wilson 
chain (calculated by the exact diagonalization method 
to highlight the discretization error) acquires a strong 
artificial temperature dependence (see Fig. 6) due to the 
discretization error. This demonstrates that even for the 
flat wide band the NRG discretization error is impor¬ 
tant to nonlocal quantities. In the flat wide band case, 
since the nonlocal Ximp is essentially a local quantity as 
Ximp = Xioc, the discretization error in Ximp could be 
largely canceled by subtracting two nonlocal quantities 
(i.e., the total susceptibilities of the system with and 
without the impurity). This validates Wilson’s defini¬ 
tion in the flat wide band. However, such a cancella¬ 
tion of the discretization error may be not always strictly 
guaranteed in arbitrarily structured bands for which the 
impurity susceptibility Ximp is a true nonlocal quantity, 
as shown in Fig. 5. This again sheds a shadow on the 
NRG calculation of the impurity susceptibility. 


D. Remarks 

It is now in the position to critically discuss the effect 
of an energy-dependent conduction band on the magnetic 
response of an Anderson impurity. The results of Sec. 
IIIB and IIIC provide a definitive answer to the prob¬ 
lem. As long as the host density of states is not gapped at 
the Fermi energy, for arbitrary energy dependence even 
though there are discontinuities or singularities in the 
band, the low-energy excitations of the impurity are al- 
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ways a Fermi liquid and paramagnetic. For such systems, 
the resultant Kondo screening of the impurity magnetic 
moment (demonstrated explicitly by the temperature de¬ 
pendence of Xioc) is qualitatively the same as in the flat 
wide band. The observation is clearly in contrast to 
the previous investigations [Ill - [l3| . This is because the 
widely used impurity susceptibility Ximp, involving the 
additional susceptibility from the conduction electrons, 
is not a good quantity to determine the impurity Kondo 
physics when the host density of states is energy depen¬ 
dent. For example, a divergent and/or negative Ximp as 
T —>■ 0 certainly invalidates the standard definition of the 
Kondo temperature Tq = limr^o 1 /4ximp W- For such 
systems, the correct magnetic properties of the impurity 
must only be drawn from the local susceptibility. 

Recently, Haul and Weichselbaum 2^ have proposed 
a new susceptibility x^'^ = Xioc + X™’’^ for the Anderson 
impurity by adding the local susceptibility xioc and the 
so-called mixed susceptibility 




(27) 


in order to achieve universal Kondo scaling for narrow 
bandwidth. This mixed susceptibility, describing the 
bath magnetization in response to a local field, is a non¬ 
local quantity and thus may be still vulnerable to the dis¬ 
cretization error. More importantly, by using the EOM 
approach, x™'’^ can be expressed in the same form of 5xc, 
i.e., as the second term in the right-hand side of Eq. (20) 
only with the magnetic field now applied on the impurity 
not on the bath. This means y™*’' can still be negative 
for some structured bands. Therefore, there is no guar¬ 
antee that in arbitrary energy-dependent bands the new 
susceptibility x’*'^ is always reliable for characterizing the 
impurity Kondo effect, even though it works well for the 
narrow-band model. 


IV. CONCLUSION 

A comparative investigation of the local and impurity 
susceptibilities for an Anderson impurity, by using the 
EDM NRG technique, has demonstrated the importance 
of the local susceptibility in characterizing the intrinsic 
magnetic properties of the impurity. Within the EDM 
algorithm, while the calculation of the impurity suscep¬ 
tibility inevitably produces severe artifacts at very low 
temperatures, the local quantity, due to its sum-rule na¬ 
ture, can be calculated very accurately in all parameter 
regimes. In particular, the accuracy of the local suscep¬ 
tibility calculated at much lower computational cost is 
at least one order higher than that of the impurity sus¬ 
ceptibility. For certain class of single-impurity Anderson 
systems in which the host density of states is arbitrar¬ 
ily energy-dependent but not gapped at the Fermi en¬ 
ergy, we have revealed that the non-Fermi-liquid and/or 


diamagnetic behaviors found in the literature based on 
the knowledge of the impurity susceptibility are spurious. 
The correct magnetic properties of the impurity in such 
systems should only be deduced from the local suscepti¬ 
bility, which suggests that the low-energy excitations of 
the impurity is always a Fermi liquid and paramagnetic. 
We hope this paper could indeed motivate more attention 
to the local susceptibility whenever the intrinsic magnetic 
response of the impurity is concerned. 
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appendix: derivation of equation (12) for the local 
susceptibility 

In this appendix, we provide the full details of deriving 
Eq. (12) within the EDM approach [13. We start from 
the spin correlation function T{t) = ([S' 2 (t), SI{0)])h ap¬ 
pearing in Eq. (11). Inserting the completeness relation 
(4) twice gives 


T{t) = Tr (^S'2/9e*"*5'2e-*"* - 

^max 

= E E 


) 


N,N'—Nmin 1,1' ,e,e‘ 
X 


{l',e';N'\S,e-^^^\l,e;N) 
(^{l,e;N\S,pei^*\l',e';N') 

(Al) 

is decomposed into 


= T^t)-tT"(t)-kT™(t), 


N,N'=N„, 


where the double sum ^ 
three contributions with N = N' (T* term), N > N' 
(T“ term), and N < N' (T™ term), respectively. The 
first contribution is 




T'{t) = E E 

N=Nraiti 1,1' ,e,e' 


{l\e'-,N\S,e-^^*\l,e;N) 


x(^(l,e-,N\S,pei^^\l',e';N) 

-{l,e-,N\pS,e^^*\l',e'-N)^ 


■ (A2) 
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The second contribution is 


T"(t) 


-^max 1 -^max 

SEE 

N'=Nmin N=N' + 1 l,l',e,e' 

x(^{l,e-,N\S,pei^*\l',e';N') 

-{l,e;N\pS,ei^*\l',e';N') 


-^max — 1 

E E 

N—Nmin k,l,e,e' 


{l,e';N\S,e-T^^*\k,e;N) 


x(^{k,e;N\S,pei^*\l,e'-,N) 

-{k,e;N\pS,ei^^\l,e';N)^ 


(A3) 


The last equality of Eq. (A3) has applied the relation (6) 
and a notation change N' ^ N, I' ^ I in the final result. 
The third contribution is 


T“(t) 


-^max 1 -^max r 

E E E 

N=N„,ir, N'=N+ll,l',e,e' 

{l',e';N'\S,e-i^*\l,e;N) 

x(^{l,e-,N\S,pei^^\l',e';N') 

-{l,e;N\pS,ei^*\l',e';N')) 


-^max — 1 

E E 




{k,e';N\S,e-i^*\l,e;N) 


x(^{l,e;N\S,pei^*\k,e';N) 

-{l,e;N\pS,ei^*\k,e';N)^ 


(A4) 


Again, the last equality of Eq. (A4) is due to the appli¬ 
cation of the relation (6). We substitute the full density 
matrix (7) into Eqs. (A2)-(A4), and then use the NRG 
approximation e; A^) « *|s, e; A^), the or¬ 

thonormality (5), and the local nature of the impurity 
spin (s, e; A^|S'z|s', e'; Af) = dee'(s; Afl-S'zis'; A^). These 
lead to 


^ max 1 


r*(t) = £ ^ 


N — Nmin k,l 




Nmax — N 




-^max 1 -^max 


-EE 


z\lk 


N^Nmin N'=N+1 k,l,e 

x{k,e;N\pN>S:,\l,e;N))-, (A6) 


-^max 1 -^max 


r™(t) = £ £ £ 

N=N,^ir,N'=N+lk,l,e ^ 


X (Z, e; A^l S'zPAf/ |fc, e;Af) 


N max — 1 


JV = JV„in k,l ^ 


N-mSix — N 




(A7) 


Here the notation [Sz]^, = (s; A^|S'z|s'; N) is in¬ 

troduced. The terms in Eqs. (A6) and (A7), 
which contain the A'th-shell density matrix pjv/ = 
Z~^^i \l,e; N'){l,e; N'\, need further calcula¬ 

tions as follows: 


£(A:,e;Af|pAr/S'z |Z,e;N) 

e 

^max 

= E E e; N\ pN, \l', e'; N") {I', e'; N"\S, \l, e; N) 

^max 

= E E e; iV| pN' \l', e'; N") {I', e'; N"\S, \l, e; N) 

N"=N ^1 V,e',e 

= £ (fc,e;iV|pAr'|fc',e';iV)(fc',e';iV|5z|;,e;iV) 

fc',e',e 

= (A8) 

k' 

Similarly, 


T* (t) 


-^max 

E E 

f=Afml„ 1,1' 




-N 



£ (/, e; iV| S.pN' \k, e;N)=J2 Pk'k [S.C ■ (A9) 


In Eqs. (A8) and (A9), the reduced density matrix 
p^f^ = ^^{k,e; N\pN'\k',e; N). We can now collect all 
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the contributions Eqs. (A5)-(A9) to obtain 

-^max — 1 -^max 

rw = 1 : 1 : 1 : 

N' — N-\-l k,k' ,l 




-^max — 1 




E St 


N — -/Vmin k,l 


^±{Er-E-)t_ME^-Er)t 


l-^max ^ 


[5.]^ [S.]fi 




+ 


Af=Afmln 1,1' 




lWmax-W 




X e 


(e--’'*.- - IS.C [SJ» }, (AlO) 


Substituting Eq. (AlO) into Eq. (11) and performing 
the Eourier transformation by using the following inte¬ 
gral representation of the Heaviside step function 


m 


1 

2TTi 



pitT 

— -dr 

T — IT] 


(All) 


with ry —>■ 0+, we obtain the retarded Green’s function of 
impurity spin 


G 5 . (e) = (e) + (e) + G^: (e), (A12) 


-^max 1 -^max 

Gs^ (^) = X! 51 X! I £ _ (]^N _ j^N\ y 

N=NrninN'=N+lk,k',l fc j + 


(A13) 


N„, 


Glie) = E E 


N=N^ 


e-{E^-En+^v 


4Afmax-Af^-lg-/3-E'' 


e-{E[^-En+^V 


iS.]u{S^]^ , 


(A14) 


Gt (^) = 


N„ 

E ^ 

N=N^ir, 1,1' 


4Afmax-JV g-/3B(f _ e-/3Bf 




e - - E^) + iri 

(A15) 


In deriving Eq. (A13), the Hermitian conditions = 
Pk'k — {^^]ki used. Substituting the above 


Lehmann representation for Gs^{s) into Eq. (10) yields 
straightforwardly Eq. (12) in the main text. 

Einally, we would like to remind the reader that the 
A^'th-shell density matrix and hence the reduced den¬ 
sity matrix p^f^ defined here are different by a factor 
^-iqWmax-Af , as compared with the original 

definitions in Ref. [17]. Nevertheless, the numerical algo¬ 
rithm for p^f^ remains the same as before. 


* Electronic address: fangtiefeng@lzu.edu.cn 
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